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I. SUMMARY 

The present research program included three major experimental studies. 

A. Spatial and Energy .Distributions of Reflected Helium Atoms 

Energj’ transfer in collisions of satellite-speed (7000 m/sec) helium 
atoms with a cleaned 6061-T6 satellite-type aluminum surface was investigated 
using the molecular-beam technique. Spatial and energy distributions of 
reflected helium atoms were measured and analysed. Energy accommodation co- 
efficients were then extracted from the measured distributions. These re- 
sults are included in a separate report (see Appendix) . 

B. Gross Accommodation Coefficient inside an Aluminum Cavity 

The gross accommodation coefficient for a satellite— speed (7000 m/sec) 

helium beam entering a 2-inch-diameter 6061-T6 aluminum spherical cavity was 

determined by measuring the exit velocity distribution of the leaving helium 

atoms using a metastable time-of-f light method. The mean velocity and energy 

reduced from the measured TOF spectrum are about 1535 m/sec and 0.054 eV 

respectively. These values correspond closely to the mean velocity and 

energy of a room-temperature effusive helium beam (l.e., ^o^en beam ” 

m/sec and E , = 0.052 eV at T = 300“K) . These results indicate 

oven beam oven 

that the 7000-ra/sec satellite-speed helium atoms entering the cavity gain 
full accommodation with the room-temperature inner surface of the sphere 
through a large number of collisions before leaving the spherical cavity. 

As will be seen in Section III, these results validate the full-accommodation 
assumption frequently used in evaluating the mean velocity of gas atoms re- 
emitted from a satellite enclosure, 

C. Arc-Heated Hydrogen Beams 

The feasibility of producing a satellite-speed atomic hydrogen beam 
by arc-heating, for use in studies of interactions of satellite-surfaces 


1 


with hydrogen atoms under laboratory conditions, was investigated. It was 
found that a stable arc-heated molecular hydrogen beam can be obtained us- 
ing Young's arc-heater, and that a partially dissociated hydrogen beam 
(i.e., a H+H2 binary beam) can be produced by using a special high- 
temperature anode. However, a pure atomic hydrogen beam was not realized. 
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II. EXPERIMENTAL APPARATUS 

The present experimental study was carried out in the UCLA Molecular- 
Beam Laboratory using the molecular beam system shovm schematically in 
Figure II-l. The satellite-speed beams were generated by an arc-heated 
supersonic-beam source developed in this laboratory. The incident beam 
was collimated by a 0,10-inch-diameter orifice placed between the collima- 
tion chamber and the detection chamber; it was characterized by a multi- 
disk velocity filter located in the collimation chamber. Since the design 
and the operation of this selector have been discussed fully in [1], only 
a brief description will be provided here. 

For a multi-disk selector, the nomina]. velocity u^ of the transmitted 
particles is determined by 


^ A 

o _ 

u u 


(II-l) 


i.e. , the transit time for a molecule with speed u^ to travel the distance 
L^ between the first and the last disk equals the time required for the 
rotor with an angular velocity OJ to rotate through the angle (f). Using a 
straight-slot approximation, one may calculate the nominal velocity of the 
transmitted molecules using 


idiere f (= ai/2ir) is the rotor frequency measured in hz and j6-5^| is a 
linear displacement related to 4> and measured in thousandths of an inch. 
(Cf. Ref. [1] for the derivation of Eq. (11-2) and the physical dimensions 
of the selector.) 


[1] S. M. Liu, SEAS Report No. UCLA-ENG-7510, UCLA, 1975. 



Figure II-1. Schematic Diagram of the Molecular Beam System. 
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Using this velocity selector, one may characterize the velocity distri- 
bution for a given incident beam by the transmitted beam signal measured by 
the mass spectrometer as a function of the rotation frequency of the selector 
at a preset displacement | 1 • The mean energy of a molecular beam with a 

large hydrodynamic speed ratio can be approximated by 

E. = 1/2 m*u^ 

X o 

where m is the mass per molecule. 

The nickle-plated detection chamber was pumped by a 10-inch oil- 
diffusion pump and a liquid-nitrogen cryopurap. The background pressure in 

_7 

the detection chamber was 10 Torr (due mostly to nitrogen and oxygen) 
under normal operating conditions. Two detection systems were constructed 
during the course of this study to accommodate the first two experiments. 

A new mechanism for measuring the complete three-dimensional spatial and 
mean energy distributions of satellite-speed helium atoms reflected from 
satellite surfaces and a new time-of-flight detector for measuring the 
velocity, distribution of helium atoms leaving a 2-inch aluminum sphere were 
built. These detection systems will be described in the following sections. 
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III. GROSS ACCOMMODATION COEFFICIENT INSIDE AN ALUMINlftl CAVITY 

An earth satellite using an on-board mass spectrometer for measuring 
the global density and composition patterns of the upper atmosphere often 
has an enclosing configuration as shown in Figure III-l. A major advantage 
of this design is that the density within the satellite is usually an order 
of magnitude higher than the ambient density, so that the sensitivity re- 
quirement of the on-board instrument can be relaxed. However, this configu- 
ration presents the problem of relating the measured density inside the 
satellite to the ambient density. It has been show that the inside density 
and the outside density are related by 


n. V ^ I T -S 

-3^ = 4*G j erf CS) [1 + + 

% V. { 2S 

X 


(III-l) 


where 


n^ = density inside the satellite 
n^ = density outside the satellite 

V = satellite velocity 

S u 

= exit velocity 

- S ^sat^\,ex^ 

V = mean random speed of ambient gas 

w,ex ^ 


G = geometrical factor (= effective entrance area/effective exit area) . 
However, in order to use Eq. III-l, one needs the value of the exit velocity 


V^. Often one assumes that particles have gained full accommodation with 
the inner surface before exiting from the satellite, so that one can deter- 
mine from the surface temperature T^ of the satellite. Possible devia- 
tions from this assumption might lead to significant errors. Hence, we 


describe here an experimental method for determining under laboratory 
conditions. 
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Figure II M. Schematic Diagram of the NASA Dual-Air-Density-Explorer Drag Satellite. 


The simplest and the most direct way to obtain the exit velocity is the 
spherical cavity method, which method simulates the actual conditions xd-thin 


the satellite. See Figure III-2. In this method, a beam with density 
and velocity enters the spherical cavity through an orifice with an effec- 
tive area The exit velocity of the atoms leaving the cavity through a 

total area exit A can be determined using one of two alternative procedures: 
o 

A. Measure the exit velocity distribution of these atoms using the 
time- 0 f “flight technique. 

B, Measure the mean exit speed from the measured steady state He 
density inside the cavity using the simplified form of Eq.(lH-l) for 


S » 1, i.e., using 
A. n, 

y. = 4 -” ‘ ~ V, 
X A n. D 

O X 


CIII-2) 


Procedure A appears more direct compared to procedure B and avoids the 
necessity for careful measurements’ of n^ and n^. Hence Procedure A was 
used in conjunction with a metastable TOP detector in this study. 

The experimental setup is sho\«i in Figures III— 2 and III— 3. A 7000- 
m/ sec incident helium beam (5/16-in. diameter), entered the 2— in. dia. 6061-T6 
aluminum spherical cavity through a 9/32-ln. dia. hole. The TOP detector 
was positioned facing a second exit hole with the same diameter. The total 
exit area is about 1% of the spherical surface area, as in the NASA Dual Air 
Density Explorer Drag Satellites. Hence, the helium atoms entering the 
spherical cavity will experience, on the average, about 100 collisions x^ith 
the inner satellite su'tface before escaping. A fraction of the helium atoms 
leaving the second exit hole xs'ere excited to metastable states (He*) by a 
pulsed electron beam. The pulse circuit for the exciter, as well as the data 
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acquisition system, are shorn in Figure III-4. Ions, a by-product of the 
excitation process, were eliminated by an electrostatic ion-deflector and 
a quadrupole mass filter. The He* TOF spectrum (i.e. , the He* signal 
measured as a function of time for He* to travel a 14-cm flight path) was 
detected by a Bendix magnetic strip electron multiplier, processed by a 
FAR waveform eductor, and recorded on a photograph and/or X-Y chart. 

Figure III-5 shows a photograph of a He* TOF spectrum and the applied 
excitation pulse displayed on a Tektronic dual-trace oscilloscope. The 
horizontal scale is 20 ps per cm. The zero— time reference is determined 
by the mid-point of the excitation pulse (i-e. , t^ == (10/2-2) ps = 3 ps). 
Figure III-6 shox^s a plot of an X-Y recorder output of a He* TOF spectrum 
expressed as a function of true flight time (i.e, , the flight time relative 
to t^)» Since the PAR v7aveform eductor has 100 memory channels, the hori- 
zontal scale is 2 ps per output channel. The mean velocity and energy of 
the helium atoms leaving the 2” aluminum sphere can be evaluated directly 
from these TOF data using the folloxiring formulas: 

a. Mean Velocity: 

i i 

b. Mean Energy: 

E - (23N^(t)(mv^^/2)/(2^N^(t)) (111-4) 

i i 

where v^ = 14 cm/(t-t^) and i corresponds to the ith channel. The detail 

calculations are given in Table XIX--1. The mean velocity and the mean energy 

so determined are about 1535 m/ sec and 0.054 eV respectively. No t^ that 

these values correspond closely to the mean velocity and energy of a raom- 

temperatiire effusive helium beam (i.e* , E == 2kT “ 0*052 eV and V = 

oven o oven 
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FILAMENT POWER SUPPLY AND 
PULSING CIRCUIT 


Figure 1 1 1-4. Block Diagram of the Electronic 



SIGNAL OUTPUT 

PROCESSOR_ 


of the Metastable Time-Of-Flight Detector. 
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Figure IH-6. Time-Of-Fllght Spectrum of Helium Atoms Leaving a 2"- Dia Aluminum Cavity. 



’ TABLE III-l 


TOP Data of Helium Atoms Leaving a 
2"-Dia. Aluminum Sphere 


T rue 

Flight 

Speed 

Signal 

Energy 

Speed 

T ime 

(ffs) 

(m/sec ) 

Height 

(eV) 

Di stri buti on 

t-to 


^i 

Ni(t) 

T.vf/2 

n.(v)% 

32 


4375 

2 

.397 

.0050 

34 


4118 

6 

.352 

.0170 

36 


3889 

3 

.313 

.0096 

38 


3684 

5 

.281 

.0177 

40 


3500 

4 

.254 

.0158 

42 


3333 

6 

.230 

.0260 

44 


3182 

8 

.210 

.0381 

46 


3043 

13 

.192 

.0676 

48 


2917 

16 

.176 

.0906 

50 


2800 

20 

.163 

.1229 

52 


2692 

25 

.150 

.1661 

54 


2593 

23 

.139 

.1648 

56 


2500 

34 

.130 

.2620 

58 


2414 

44 

.121 

.3637 

60 


2333 

43 

.113 

.3804 

62 


2258 

51 

.106 

.4818 

64 


2188 

59 

.099 

.5939 

66 


2121 

62 

.093 

.6637 

68 


2059 

65 

.088 

.7386 

70 


2000 

70 

.083 

.8429 

72 


1944 

74 

.078 

.9427 

74 


1892 

81 

.074 

1.090 

76 


1842 

82 

.070 

1.164 

78 


1795 

84 

.067 

1.256 

80 


1750 

86 

.063 

1.352 

82 . 


1707 

86 

.060 

1.421 

84 


1666 

85 

.058 

1.474 

86 


1628 

88 

.055 

1.599 

88 


1591 

86 

.052 

1.637 

90 


1556 

81 

.050 

1.612 

92 


1522 

87 

.048 

1.810 

94 


1489 

86 

.046 

1.867 

96 


1458 

79 

.044 

1.789 

98 


1429 

82 

.042 

1.935 

100 


1400 

81 

.041 

1.991 

102 


1373 

77 

.039 

1.969 

104 


1346 

78 

.038 

2.073 

106 


1321 

73 

.036 

2.016 

108 


1296 

65 

.035 

1.863 

110 


1273 

65 

.034 

1.933 
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TABLE III-l (CONT'D) 


t-t 

1 

0 


Ni(t) 

mv|/2 

(v 

112 

1250 


60 

.032 

1.850 

114 

1228 


65 

.031 

2.076 

116 

1207 


64 

.030 

2.116 

118 

1186 


53 

.029 

1.814 

120 

1167 


51 

.028 

1.805 

122 

1148 


51 

.027 

1.865 

124 

1129 


46 

.026 

1.738 

126 

1111 


35 

.026 

1.366 

128 

1094 


40 

.025 

1.611 

130 

1077 


39 

.024 

1.620 

132 

1061 


39 

.023 

1.670 

134 

1045 


39 

.023 

1.721 

136 

1029 


38 

.022 

1.727 

138 

1014 


34 

.021 

1.591 

140 

1000 


32 

.021 

1.541 

142 

986 


25 

.020 

1.239 

144 

972 


28 

.020 

1.427 

146 

959 


25 

.019 

1.310 

148 

946 


25 

.019 

1.346 

150 

933 


20 

.018 

1.106 

152 

921 


21 

.018 

1.192 

154 

909 


24 

.017 

1.399 

156 

897 


22 

.017 

1.316 

158 

886 


23 

.016 

1.411 

160 

875 


16 

.016 

1.007 

162 

864 


13 

.015 

.8384 

164 

854 


10 

.015 

.6610 

166 

843 


16 

.015 

1.083 

168 

833 


14 

.014 

.9710 

170 

824 


10 

.014 

.7102 

172 

814 


16 

.014 

1.163 

174 

805 


16 

.013 

1.190 

176 

795 


12 

.013 

.9135 

178 

787 


13 

.013 

1.012 

180 

778 


8 

.013 

.6370 

182 

769 


11 

.012 

.8954 

184 

761 


6 

.012 

.4992 

186 

753 


10 

.012 

.8502 

188 

745 


8 

.012 

.6949 

190 

737 


A 

.011 

.3549 


(t) = 3247 





iNi 

(t).Vi/£N.(t) 

= 1636 

m/sec 




(t) . (l/2)mv^ / 

ZN,(t) 

= 0.054 

(eV) 
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1550 m/sec at T “ 300®K). Hence it is concluded that the 7000-m/sec 
oven 

incident helium atoms entering the aluminum spherical cavity are fully 
accommodated with the room-temperature inner surface through a large numher 
of collisions (about 100) before escaping from the cavity. 

The velocity distribution can be obtained from the TOF spectrum us- 
ing 

n^(v) = N^(t)*(t-t^)^/(]2N^(b)'Ct-t^)^) 

' ' i 

The results are shorn in the last column of Table III-l and in Figure 


III-7. 



IV. ARC-HEATED HYDROGEN BEAMS 


The feasibility of producing a satellite-speed hydrogen beam using arc- 
heating was investigated. More specifically, the design parameters, start- 
ing procedures, and operating parameters for obtaining a stable (and repro- 
ducible) satellite-speed hydrogen beam which would facilitate studies of 
interactions of hydrogen atoms with satellite- type surfaces were studied. 

The arc-heater used in the scattering experiments to produce satellite- 
speed helium beams was used, with minor modifications, for this study. In 
particular, the nozzle diameter was increased from 0.050" to 0.100" to prevent 
melting of the copper anode due to excessive heat load. It was found that 
the most critical operating parameter is the hydrogen flow rate; the hydrogen 
arc is stable only within a narrow range of flow rate. If the flow rate is 
too high, the nozzle is eroded due to over heating. If the flow rate is too 
low, the arc becomes internal and damages internal components of the arc- 
heater. 

It was found that a stable arc-heated hydrogen beam can be obtained 
using the following design parameters, starting procedure, and operating 
conditions : 

A. Design Parameters 

1. Nozzle: a threaded copper anode (replaceable) with a 0.100"-diameter 

orifice. Such a nozzle can be used repeatedly. 

2. Anode- Cathode Gap; approx. 0.125" 

3. Input Current to tlagnetic Coil: approximately 10-12 amps. 

Field Strength: 600 Gauss 

B. Starting Frocedure and Operation Conditions 

1. Light the argon arc -sri-th high argon flow rate. The power input is 
approximately 120 amps/15 volts. 
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2. Adjust P to 50 torrs by lowering the argon flow. The input 

stg 

power might shift to 120 amps/20 volts. 

3. Add hydrogen slowly until hydrogen flow rate reaches 10 mole/ 
sec (22 std cc/sec). The stagnation pressure decreases first then 
increases beyond 50 torrs. It is necessary to use a good needle 
valve to control hydrogen flow. 

4. Turn off argon flow slowly; keep input current at 130 amps. The 

input voltage will increase to approximately 30 volts. ^^stg^hot ~ 

35-37 torrs. (P = 8.5 torrs.) 
cold 

5. The input voltage will drop to 26-27 volts after a few minutes of 
running. The input current can be varied between 90 to 150 amps. 

The hydrogen flow rate can also be lowered to obtain various 
operating points. 

Hydrogen beams produced using these techniques were characterized by 
the multi-disk velocity selector described in Section II. Typical measured 
velocity distributions are shown in Figure IV-1. The results show an excel- 
lent linear correlation between displacement |S-6^| and frequency f^ corres- 
ponding to the peak of the velocity distribution, as predicted by Eq. II-2. 
The approximate nominal velocities and corresponding beam energies are also 
shown in Figure IV-1. 

Way et al. [2] reported some success in producing an atomic hydrogen 

15 2 

becm with a flux on the order of 10 atoms/cra -sec at a distance of 1 meter 
from the source. They also used Young's arc-heated beam-source design with 
minor modifications. Although they reported that energies of their atomic 
hydrogen beams ranged from 1.3 to 5.5 eV, the uncertainties in their measure- 
ments and the results for the molecular hydrogen beam obtained in this 


[2] K. R. Way, et al., Int. Conf. Phys. Elect. Atomic Collisions , Seattle, 
Washington, July 1975. 
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laboratory lead us to believe that the energy of the atomic hydrogen pro- 
duced by this sourse is less than 2 eV. Winklemann [3] also produced a 
H+H 2 beam (with 30% dissociation) using Young’s arc-heated beam-source 
design with a special W-Cu anode and a 0 . 050”-diameter orifice. 

From the results obtained in this study and information provided by 
other investigators, the following conclusions regarding the production of 
arc-heated hydrogen beams are drawn: (1) a molecular hydrogen beam can 

be obtained using the design and procedures described at the beginning of 
this section; (2) a partially dissociated hydrogen beam (i.e. a 
nary beam) can be produced by using a special high-temperature anode which 
can end^^rc a higher heat load and therefore allows the source to operate 
with a smaller nozzle orifice and a lower hydrogen flow rate; and (3) a pure 
atomic hydrogen beam would be difficult to obtain using arc heating. 


[3] Private communication. He is at the Max-Planck Institut fUr 
StrHmungsforschung, GSttingen. 
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SPATIAL AHD ENERGY DISTRIBUTIONS OF 
REFLECTED HELIUM ATOMS 



FOREWORD 


The researches described here were supported mainly by the National 
Aeronautics and Space Administration (under Grant NGR 05-007-416) and by 
the UCLA School of Engineering and Applied Science* These studies were 
part of a continuing program of researches In gas-surface interactions* 


ABSTRACT 


Energy transfer in collisions of satellite-speed (7000 m/sec) helium 
atoms with a cleaned 6061-T6 satellite-type aluminum surface was investi- 
gated using the molecular-beam technique. The amount of energy transferred 
was determined from the measured energy of the molecular-beam and the 
measured spatial and energy distributions of the reflected atoms. 

Spatial distributions of helium atoms scattered from a 6061-T6 aluminum 
surface were measured again in this study, and show features similar to 
those presented in report UCLA-EKG-7546 [1] . The scattering pattern exhibits 
a prominent backscattering , probably due to the gross surface rougnesss and/ 
or the relative lattice softness of the aluminum surface. 

Energy distrlbutiphs of reflected helium atoms from the same surface 
were measured for six different incidence angles. For each incidence angle, 
distributions x^ere measured at -approximately sixty scattering positions. 

At a given scattering position, the energy spectra of the reflected helium 
atoms and the background gas were obtained using the retarding-field energy 
analyzer. The mean reflected-beam energy and the differential energy accom- 
modation coefficient ((A.C. ) were then extracted from these 

spectra using a least-square fitting program. The measured (A.C.)„C0.,0 ,4*)s 

b X r 

shoxT some fluctuations and a weak dependence oh scattering angle, i.e., the 
accommodation decreases sloxofly as the scattering direction shifts toward 
the surface tangent. 

The overall energy accommodation ioefficient for a beam x^ith a given 
incidence angle was then evaluated using the measured spatial density dis- 
tributions and the mean reflected-beara energy distributions. Results shoxf 
that the mean accommodation coefficient varies betxjeen 50% and 65%, depen- 
dent on the incidence angle. 
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LIST OF SYtmOLS 


(A.C. differential energy acconuuodation coefficient at scatter- 
ing position C0^,(})) for a beam with 0^ incidence angle 

(A. 00^(0.) overall energy accommodation coefficient for a beam with 

0^ incidence angle 

incident -beam energy 

E C0.J0 energy of scattered helium atoms at scattering position 

^ ^ ^ (0^j^) a beam with 0^ incidence angle 


E 


r 


energy of scattered helium atoms evaluated from its dif- 
ferential energy spectrum 


f(E) 

n^ce^.e^.W 


mean reference energy of thermal background gas evaluated 
from its differential energy spectrum 

differential energy distribution 

normalized spatial density-distribution function 

surface normal 


surface tangent 

0. incidence angle of helium beam measured from surface 

normal 


0^ in-plane scattering angle measured from surface normal 

(j) out-of-plane scattering angle measured from the plane of 

incidence 


ff 


standard deviation of the true reflected-beam energy 
spectrum from the least-square fitted function 
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CIIAPTER I 
INTRODUCTION ' 

) Basic knowledge concerning energy and momentum transfer between earth 

satellites and upper-atmospheric gases is essential for understanding the 
drag experienced by earth satellites (therefore for estimating the life- 
time of an earth satellite and/or extracting the mean upper-atmosphere 
density from satellite drag data). For example, in predicting the aero- 
dynamic drag of a satellite, one uses frequently a model in which the 
thermal accommodation between the ambient gas and the satellite surface 
is complete and the scattering distribution of reflected molecules follows 
the cosine law. However, possible deviations from this model might yield 
greatly different results. 

These energy and momentum transfers can be investigated experimentally 
in the laboratory using an ultra-high vacuum system and the molecular-beam 
technique. The desired information can be extracted from the change in the 
beam properties during the surface collision if the states of both the inci- 
dent anfl the scattered beam (spatial distribution and speed distribution) 
can be determined. Spatial distributions of satellite-speed helium beams 
scattered from satellite surfaces were obtained previously and suramariaed 
in report UGLA-ENG-7546. This report presents measured energy distributions 
of helium atoms reflected from 6061-T6 aluminum surfaces. 

In Chapter II, the experimental apparatus and procedures are described 
briefly. Emphasis is given on the design and the operating procedure for 
the retarding-field energy analyzer. Experimental results are given and dis- 
cussed in Chapter III. A least-square curve-fitting computer program is 
given in an Appendix. 
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CHAPTER II 


EXPERItlENTAL APPARATUS AND PROCEDURES 

The present experimental study was carried out in the UCLA Molecular- 
Beam Laboratory using the molecular-beam system shovm schematically in 
Figure II-l. Since it has been described in detail elsewhere [1,2], only 
a brief description tfill be given here. 

The satellite-speed (7000 m/sec) helium beams were generated using an 
arc-heated supersonic beam source developed by Young [3]. The incident 
beam was collimated by an orifice of 0.10-inch diameter placed between 
the colliraation chamber and the detection chamber. The beam was character- 
ized by a multi-disk velocity selector located vn the collimation chamber. 

A new detection system was constructed during the course of this study 
for facilitating measurements of the complete three-dimensional density and 
mean-energy distributions of satellite-speed helium atoms reflected from 
satellite surfaces. Cf. Figure II-2. This net/ system includes (1) a target 
positioning mechanism, (2) a detector rotating mechanism and (3) a mass 
spectrometer atid/or a retarding-field energy analyzer. Descriptions of the 
first two mechanisms were given in the first report of this study (cf, ref. 1) . 
The design and the operating procedure for the retarding-field energy ana- 
lyzer t/ill be given here. 

The retarding-field energy analyzer is shown in Figures II-3 and II-4. 

An electron- impact ionizer, mounted 0.5-inch from the target surface on the 
entrance plate of the analyzer, was used to ionize a fraction of the beam 
species (also of the residual background) , The retarding-field section of 
the analyzer is made of seven thin stainless-steel washer-shaped discs 
placed in a stainless-steel can. The inlet plate is followed by three 
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Figure IM, 



Schematic Diagram of the Molecular Beam System. 






















Figure ||.3 Photograph of the Retarding Field Energy Analy 
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Figure II4» Schematic Diagram of the Retarding-Field Energy Analyzer. 



focusing plates, a retarding plate and two more focusing plates. The 
potentials of all plates except that of the inlet plate are floated rela- 
tive to the retarding potential so that ions of different energies will 
experience the same focusing effects within the analyzer. Thus undesir- 
able effects of the energy-dependence of the transmission efficiency 
are minimized. Typical plate potentials also are given in Figure II-4. 

The ions that have passed through the retarding-field region were 
filtered by a 2-inch quadrupole mass filter to eliminate the noise from 
the ionized background gases. The filtered ions were then detected by a 
pulse-counting particle detector. The energy spectrum of the reflected 
atoms at a given scattering position was obtained by measuring the 
reflected-beam density as a function of the retarding potential. The 
measured spectrum was processed by a KS513 signal averager and recorded 
on IBM cards. A block diagram of the electronic system is shown in Figure 
II-5. 

Although the electron-impact ionization does not change the kinetic 
energy of a helium atom (since the translational energy transfer between 
the ionizing electron and the atom is negligible due to the large ratio of 
their nasses) , it was found that space-charge effects of ionizing electrons 
in the ionization region and/or surface-charge effects on the anode cage 
did introduce a systematic shift of the entire energy spectrum toward lower 
energies (i.e. , the positive ions were produced in a region of negative 
potential with respect to ground). To reduce this shifty a small emission 
current (y 50 yA) was used in the ionizer. Also, a positive potential (8 
volts relative to ground) was applied to the anode cage in order to counter 
shift the energy spectrum toward higher energies. Then, since the potential 
of the ionization region was no longer at ground level, it was necessary to 
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Figure II-5. Block Diagram of the Electronic System of the Retarding-Field Energy Analyzer. 
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ascertain independently a reference point in the energy spectrum. This 
reference paint was provided by the thermal energy spectrum of the back- 
ground gas, x^hich has a mean thermal energy of 0.05 eV (i.e,, the mean 
thermal energy at 296°K) . 

Since the background gas of the beam species also contributed to the 
measured spectrum, it was necessary to subtract this contribution in order 
to obtain the reflected-beam energy distributions. This subtraction xjas 
facilitated by measuring tX70 spectra (one for the reflected beam plus back- 
ground and one for the background alone) under the same operating conditions 
Both spectra were then processed using a computer program; the reflected- 
beam energy spectrum was obtained by subtracting the background spectrum 
from the overall beam-energy spectrum. Both the background spectrum and 
the reflected-beam spectrum were least-square fitted using a high-order 
Chebyshev polynomial function. The differential energy distributions f(E) 
x^ere obtained by simple differentiation of the fitted functions. The mean 
reflected-beam energy at a given scattering position xras evaluated from 


+ 0.05 (eV) (II-l) 


x^here 


E(6.,0 , 
X r‘ 


4 >) 




f(E)*dE 


(11-2) 


and 0.05 eV is the thermal energy of the background gas at 296 “K. The dif- 
ferential energy accommodation coefficient at a given scattering position 
x«7as obtained using 


E,~E (8.,0 ,(j)) 
1 r X r 


(II-3) 


[A.C. ]„(0 . , 0 ,<j)) — 


where is the incideat-heam energy. The computer program and its input 
parameters for handling the described data reduction are given in the 
Appendix. The overall energy accommodation for a given incidence angle 
was then evaluated by 




9 d) 
r ^ 


(n-4) 


where njj^(0^,0^,$) is the normalized spatial density-distribution function 
of reflected helium atoms* 

As indicated, spatial distributions of satellite-speed helium beams 
scattered from four different satellite surfaces were obtained in the first 
phase of this study. Experimental procedures and results are included in 
report UCLA-ENG-7446 [1] • However, spatial distributions at some angles in 
the backscattering region were not measured at that time due to the con- 
straint on the detector path as indicated in Figure II-6-(a). This problem 
was solved later by rotating the surface counterclockx<rise beyond the normal 
incidence angle while retaining the previous detector path as shoxm in 
Figure II-6-(b) . Spatial and energy distributions were measured in the 
present study using these complementary configurations. 
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DETECTOR PATH 

{5N THE PLANE OF INCIDENCE) 


CONFIGURATION (a) 



CONFIGURATION (b) 

Rgure ri-6. CompieiTientary Beam-Surface Confiaurations of the Scattering System. 


40 


CHAPTER III 


RESULTS AND DISCUSSIONS 

Spatial distributions of a satellite— speed (7000 m/sec) helium beam 
scattered from a cleaned 6061-T6 aluminum satellite surface for six dif- 
ferent incidence angles (0“, 15“* 30°, 45°, 60° and 75° from the surface 
normal) are shovm in Figures III-l to III-6. The center of the polar dia- 
gram corresponds to the point of impingefiient . The incident beam impinges 
on the test surface (which coincides with the surface of the page) from 
the bottom of the diagram with the given incidence angle measured from the 
surface normal. The upper (0^ > 0) and lower (0^ < 0) halves of the dia- 
gram represent the forward-scattering and bactarard- scattering regions 
respectively. The dashed lines at constant value of 0^ indicate detector 
paths (i.e. , from <j) = 0° to = 90°). denotes the out-of-plane scatter- 
ing angle and <}> = 0° represents the plane of incidence. These results show 
diffusive scattering patterns and exhibit trends similar to those previously 
reported [1] . As indicated before, the most interesting feature on these 
scattering patterns is the prominent backscattering of the incident helium 
atoms (i.e., a large fraction of the incident atoms are scattered back in 
the vicinity of the incident beam) , particularly as the incidence angle 
increases toward the surface tangent (i.e., for large values of 0^). This 
large fraction of backscattering could be due to the gross surface roughness 
and/or the relative lattice softness of the aluminum satellite surfaces. 
Smith [4] observed a large increase in backscattering Intensity for increas- 
ing surface roughness in his computer simulation of gas molecule reflections 
from rough surfaces. This backscattering could result in relatively high 
drag coefficients for such satellite surfaces. The spatial-distribution 




,0 SYMBOL DENSITY 

i (ARBITRARY UNITS) 



• 1 

0.5 


A ' 

1 

1.0 

-45“ 

O i 

1.5 

□ 

2.0 


o . 

2.5 

-60“ 

A 

3.0 


V 

3.5 

-75“ 

B 

4.0 


Figure 111-2. Polar Plot of Scattered -Beam Density Distribution for 7000 m/sec i 
Helium Beam Scattered from Cleaned 6061 -T6 Aluminum Plate at ‘ 
15® Incidence Angle. i 






SYMBOL 


DENSITY 




(ARBITRARY UNITS) 


• 

0.5 

-45° 

A 

1.0 


O 

1.5 


□ 

2.0 

-60° 

o 

2.5 

0 

▲ 

3.0 

-75 

V 

3.5 


■ 

4.0 


Figure !II4. Polar Plot of Scattered -Beam Density Distribution for 7000 m/sec 
Helium Beam Scattered from Cleaned 6061-T6 Aluminum Plate at 
45® Incidence Angle. 







Figure II 1-5. Polar Plot of Scattered-Beam Density Distribution for 7000 m/sec 
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measurements shorn here and the energy-distribution measurements to be pre- 
sented next provide the data required for estimating the overall energy 
accommodation for this beam-surface combination. 

Measurements of energy distributions of satellite-speed helium atoms 
scattered from a cleaned 6061-T6 aluminum satellite surface were made for 
six different incidence angles (6^ = 0“, 15®, 30®, 45“ , 60“ and 75“ from 
the surface normal). For each incidence angle, distributions were measured 
at approximately sixty scattering positions. These scattering positions 
included eleven in-plane scattering angles (6^ = ±75“, ±60“, ±45“, ±30“, 

±15“ and 0“) and six out-of-plane scattering angles ((|) = 0“ , 15“ , 30“ , 45“ , 
60“ and 75“ ). Typical energy spectra obtained at a given scattering posi- 
tion are shown in Figures III-7 and III-8. Curve A of Figure III-7 repre- 
sents the energy spectrum of a reflected helium beam superimposed on an 
energy spectrum of the background helium gas. Curve B of Figure III-7 
represents the (thermal) energy of the background helium gas (mostly due 
to beam load) . The reflected beam energy spectrum is then the difference 
of these two spectra (i.e., A-B). Figure III-8 shovis the normalized energy 
spectra of the thermal background and the reflected helium atoms (Curves 3 
and 1), their least-square fitted curves (Curves 4 and 2) and the correspond- 
ing differential energy distributions (Curves G and A) obtained using the 
computer program shoxm in Appendix A. The differential energy accommodation 
coefficient was obtained using Equations (II-l) and (II-3) . Results for 
(A.C.)-p(0 . ,0 ,(fi) obtained at all possible scattering angles are given in 
Tables III-l to III-6. Measurements were not possible within a solid angle 
around the incident beam (due to interference between the detector and the 
incident beam at these scattering positions) and for some glancing scattering 
angles (due to weak signal-to-noise ratios) . These tables also include 
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Table HI-1. The Differential Energy Acconnmodation Coefficients and the Normalized Spatial 

Density Distribution for 7000 m/sec Helium Beam Scattered From Cleaned 6061 -T6 
Aluminum Plate at 0^ Incidence Angle 
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(a) The Differential Acconnmodation Coefficient {%) 

(b) The Normalized Spatial Density {%) 

(c) Standard Deviaiton {%) 














































NOTE: 

(a) The differential Accommodation Coefficient (%} 

(b) The Normalized Spatial Density (%) 

(c) Standard Deviation [%) 
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Table III-4. The Differential Energy Accommodation Coefficients and the Normalized Spatial 

Density Distribution for 7000 m/sec Helium Beam Scattered from Cleaned 6061 -T6 
Aluminum Plate at 45° Incidence Angle 



IMOTE: ; 

(a) The Differential Accommodation Coefficient (%) 

(b) The Normalized Spatial Density {%) 

(c) Standard Deviaition {%) 
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NOTE: 


(a) The Differentia! Accommodation Coefficient {%) 
{b) The Normalized Spatial Density {%) 

(c) Standard Deviation {%] 
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NOTE: 
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(c) Standard Deviation {%) 
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standard deviations (cr) of the reflected-beam energy-spectrum data from the 
least-square fitted curves and the normalized spatial-distribution function 
of the reflected helium atoms obtained from the measured spatial distribu- 
tions shown in Figures III-l to III-6. The overall energy accommodation 
coefficients at a given incidence angle was then evaluated using Equations 
(II-4) and the data given in these tables. The results are shown in Figure 
III-9. 

The differential accommodations obtained show some fluctuations, due 
perhaps to the weak signal-to -noise ratio which results from the relatively 
diffusive scattering from the satellite-type aluminum surface. The results 
also indicate a weak dependence of accommodation on scattering angle, i.e. , 
the (A.C.)g(0^,9^,(j)) decreases as the scattering direction shifts toward 
the surface tangent. 

The overall accommodation coefficient is slightly higher for a glancing 
incident beam than for a normal incident beam. The value varies between 


50% and 65% for this beam-surface combination. 


OVERALL ENERGY ACCOMMODATION COEFFICIENT (%) 
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APPENDIX 


COMPUTER PROGRAM FOR LEAST-SQUARE FITTING 


A. PROGRAM : 


WAIN program: input; call; output 

CIMENSICN NC(SI2), 

I YC512) ,YBM(512) ,YaG(512) f 

1 DYLS(£T2>, 

" 1 EM2C2 J ,SIGMA2(2) , ' 

1 GF(2£1,101) ,GV(251) 

COMMON /AAA/X(512)*YNI£12) » YLS ( 5 1 2 ) , DY t 5 1 2 ) » N » £M . S IG MA . KM 
DATA BLANK/IH /»D0T/1H*/ 

DATA DM/IH 1/ ,BML/lH2/',CEM/lriS/ 

DATA aG/lH3/»BGL/lHA/»CBG/lHG/ 

i CONTINUE 

READ (5,600 J NSET , KMCN , KCA Y, KYEAR » THETA I , THETAR , PH I 
IF (NSET.LT. 1) GO TO 1000 

IvRITE (6,900 ) KMON,KCAY,KYEAR.THETAI ,ThETAR,PHI 
IF (NSET,EQ.2) GO TO 5 
IF (NSET.GT,2) NSET=2 

READ (5,810) EIiM,NCVl,Vl,NCV2,V2,NI ,XLL,NF,XHL, 

1 INTI, INTF,KPG,KSTEP,NPV;,NSTEP 

K=NF-NH-1 

DV={ V2-VU/C NCV2-NCV 1 + 1 ) 

VI=V1 +0V»{ NI-NCVl ) 

VF=V1+DV*( NF-NCVl ) 

5 CONTINUE 

IF (KPG.LT. 13 GC TO 35 
CO TO 1=1,251 
GVC I ) = 0 
DO 10 J=l,101 
GF( I, J)=t!LANK 
10 CONTINUE 

DO 20 1 = 1 ,251 ,5 
CO 20 J=l,101 
GF(I,J)=DCT 
20 CONTINUE 

CD 30 1=1,251 

DO 30 J= 1,1 01, 1C 
GF(I,J)=DOT 
■ 30 CGNTINUE “ ' ' 

55 CONTINUE 

CO 40 1=1,32 
K=16*( 1-1 ) 

40 READ { 5,830) TYBG(K+J) ,J = 1',16) 

IF (NSET,LTa2) GO TO £5 
DO 50 1=1,32 
K=16*(I-n 

READ "(5,830) { YEM ( K + J ) , J =1 , 1 6 ) 

DO 50 J=l,16 

50 YdM (K+ J )=YBM ( K+J )“YBG (K+J) 

55 CONTINUE 

CO 500 KI=1 ,N£ET 
DO 60 I=M * NF 
K = I -M+1 
NC C K) =I 
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non 


X(K)=VI+DV/'Sf(K-l ) 

GO TO (61 ,62) ,KI 
ei Y{K)=VBG(I> 

GO TO 60 
62 V{K)=Y8M(I) 

60 CONTINUE 
' SYMAX=0«0 
SYMIN=0«0 
DO 70 1=1,50 
SYM1N=SYKIN+Y (N-I ) 

70 EYMAlf^SYMAX+Y^ n 
YMIN=SYMIN/50,0 
YMAX=SYWAX/50 #0 
Ot 80 K = 1,N 

"■ ' YN ( !<■)■= ( Y { XT- YMl NT? ( YM A X- YM I N ) 

IF {YN(K).GE,0,0) GO TO 75 
YN( K)=0,0 

75 IF (YNCK)»LE.1.0 ) GO TO 80 

YNrio = i;o ■ ... 

so CONTINUE 

CA LL L EAST-SQLAFE CURVE FITTING 

CALL LSFCHV (N,XLL, XHL ) 

GO TO (110, 120) ,K I 
110 WRITE ( ,910) 

• - “ GO-TD'TSO" 

120 WRITE (6,515) 

130 CONTINUE 

IF (NPW.LT.l) GO TO 200 
WRITE (6,920) 

CO 190 K=1,N,NSTEP 

WRITE (6,910) NC ( K ) , X (K ) ,Y (X) , YN (X ) , YL3 (K ) . DY (.< ) 
190 CONTINUE 
200 CONTINUE 
EiV=0»0 
SIGMA=0,0 

CALL EMEAN( INTI , INTF ) 

EM2(KI)=EM 
CALL DEVIA 
SIGMA2(KI)=SIGy,A 
SIGNAP=SIGMA*1 00,0 
WRITE" { 6, 970 ) EM,SICyAF 
IF (KPG,LT,1) GO TC 400 
KIN=1 

CO 290 K=1,N.KSTEF 

KIN=KTN+l 

GV(KIN)=X(K) 

JYN=1004A3S (YN(K ) ) +l 

IF (JYN.LT.lCl) GO TC 291 

JYN=^t01 

291 CONTINUE 

IF (YLS(K).6£,0»0) GC TO 294. 

JYL3=1 

GO TO 292 
294 CONTINUE 

JYLS = 1 004 (AISTYLS(K) ) ) + l 
IF ( JYLS.LT.lOl ) GO TO 292 
JYL3=101 

292 CONTINUE 
JDY=1004AQS(OY(K) )+l 
IF (JDV.LT.lOl) GO TC 293 
JDY=101 

293 CONTINUE 
GO TO (300,310}.KI 

300 GF(KIN, JYN ) =oG 

GF(KIN, JYLS)=6GL 
GF(KIN, JDY)--OBG 
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I 
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f 
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c 

c 

c 


GO TO 235 

310 GF(KIN,JYN)=8M 

GF( KIN* JYLS)=BML 
GF( KINS JCY }=CBM 
235 continue 
230 CONTINUE 

IF (NSET.GT.l .AND.KI .LT«2) GO TO 400 
WRITE C6»940) 

" GV(1)=:GVT2)~(GV(3)-GV(2)) 

DO 360 1=1, KIN 

WRITE (6,350) GV C I ) , ( GF ( I . J ) » J=1 , 1 01 ) 
360 CONTINUE 

400 'CONTINUE 

500 CONTINUE 


CALCULATE_jrRJJE 


ACC=0*0 

IF (NSET,LTb 2) GO TO 500 
ETRUE=EM2( 2)-EM2{ 1 )+0,05 
~ '"A‘C'C=(Er-ETRUE1/EI*iOO ' 

600 CONTINUE 

WRITE (6,975) 

WRITE (6,980) ETRUE,ACC 
~ 'WRITE (6,373') 


GO TO 

aCC FORMAT 
810 FORMAT 
8 3CTF0RMAT 
910 FORMAT 
900 FORMAT 
1 

915 FORMAT 
920 FORMAT 
1 

93C FORMAT 
940 FORMAT 
1 
2 

950 FORMAT 


( 41 10,3F10,3 ) 

(F3,2, I3,4< I3,F7,2) ,615) 

(leFS-.o) ' — ' 

(12X,23Hy=5s SET i: BACKGROUND *«/) 

(///2X .7FDATE : • 12 , 1 H/ , I 2 , IH/ , 1 2 , 5 X, 

26HANGLES(ThETAl/ThETAR/PHI ) : , F6 . 2 , 1 H/ , F6, 2 « 1H/,F6<,2/) 
■(12X,17H:f:* SET 2i BEAM '^*/) 

(3X,7FCH, NO, , 12X,3HV-R,1 1X,5HI-SIG, 13X,SHI-NJ0R, 
10X,7HLSI-N0P. ,ax,9HD{LSI)/0V) 

(5X,IS,1CX,F6*3,1CX,F6. 1,3( lOX. F7.4) ) 
(//8X,3HV/R,10X,7FI : I (BM) i3X, 9H2:LSr(aM), 3X. 
loHB ;D{ LSI ) /D V-( 8M) ,3 X , 7H 3 : I C UG ) , 3X , 9 F4 :LS 1 ( BG ) , 

3X, 16HC(LS I )/DV~( EG )//) 

(2X ,F1C.S,3X, 101 A1 ) 


370 FORMAT ( /5X , 12HMtA N-ENEKGY : , F7 .4 ,3 ri EV » lOX , 19riSTANDARD DEVIATION:, 
I 2X,F5.2, IHX/) 

37S FCRiMAT = 

1 of* ** :{, -je^c# j(: ^ ^:s: // ) 

9dO FORMAT (/lOX, 22HTRUE MEAN-OEAM EN ERGY : , 2X ,F7 , 4 , 3X , 2HE V , 1 0 X , 

1 26HACCQMMaDAT I CN CCEFFI CI£NT= , F7, 4, IH 5 /^ ) 

C 

1000 CONTINUE 
.STOP 
END 


C 

C 


c 


c 


SUBROUTINE LSFCHV { M , XLL , XHL) 

COMMON ^/AAA/XlSlil^YCSl 2), YLSt 512) ,DY( 512) , N , EM . SI GMA , KM 
COMMON /SIM/A( 100 ) ,R C 10 } . . 


IMTIALI2ATICN 
DO 2 I =1 iMl 
2 R( I )=0 »0 
N2=M1 =M1 


DO 4 1=1, M2 
4 A( I ) = 0 .0 

X0=2«0/(X(N)-X( 1) ) 

XO=(X( N)+X(l ) )/CX(N) -X( 1 ) ) 
CALCULATE ChEBYSHEV FUNCTION AT 
T( 1 )=1 .0 
DO 30 K=1,N 
TX=X04X(K1-X0 
T(2)=TX 
DO 10 I =2, Ml 

10 T( I ) =2 •0’!'TX*T{ I-l ) -T { r-2 ) 
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Dd 20 1=1. Ml 

R( I )=R H )+T( I )«Y(K ) 

DD 20 J=1.M1 

IJ = ( I-l J=!=M1 + J 

20' A£ ij)=A'( rji+T( n^T( j ) ■ ■ 

30 CONTINUE 

EPS=1 . OE-60 
IER=0 

call SIMC(M1 ,1 .EPS, IER) 

calculate ChEBYSHEV PCLYNCMIAL AT TX USINij R(I) 
CALCULATE OY FORM CHEEYShEV POt_yN(JMIAL 
DT( 1 ) =0.0 

"CT { 2 ) =r'IO ■ 

DQ 30 K=1,N 
YLSK=0,0 
DYK=0 .0 

Tx=xo»'xrto-xo — 

T{ 2 )=TX 
CD 3S 1=3, Ml 

35 T( I }=2,0=fTX«T(I-1 }-T C I~2 3 

D'a"4CT"l = r7Ml ■ 

40 YLSt<=YLSK+T( I )*R(I ) 

50 YLS{K)=YLSK 
DO 60 1=3. Ml 

~ 5-0'“DT ( I ) = 2T0-^Trn'-T) +TX 4 CT £ I-ll ) -DT £ 1-2 ) 

DO 70 1 = 2, Ml 

70 CYK=DYK+Rt I ) =iDT £ n 

IF £DYK.LE.0.03 GO TC 80 
DYK = 0* CCCOO'Ol 
aO CYCKJ=DYK 
DYM=0. 0 
CO 250 K=1,N 

IF £X£K}oLT.XLL*OR-,X£K).GT.XHL) GO TO 2oS 
240 IF £DY£K)oGT»0.0 ) GO TO 250 

■~IF £ DYM.GT.ABS£DY( K) ) ) GC TO 250 
CYM= ABSIDY (K ) ) 

KW = K 

250 CONTINUE 
■KRS=0 " ■ 

K=0 

230 K=K'1-1 

KR=KM-K+1 

TF £KR.LT.l) GO TO 2R0 
IF £KRS.GT,1) GO TO 2E5 
CY £KR)=DY£KR)/OYM 
IF £DY£ KR) .LT.O .0 ) GC TO 230 
KR S “2 

2£5 CY £ KR> =-0.0C OOOlXDYM 
GO TO 280 
290 K=0 

KFS=0 ' 

292 K=K+1 

kf=kmfk 

IF (KF.GT.N) GO TO 300 

IF (KFS.GT.l) GO TO 254 

DY £ KF)=DY£ KF >/OYM 
IF £ DY £KF> .LT.0.0 1 GC TC 292 
KFS=2 

■294”DY £ KF )=-C.0C OCOT/D'VM 
GO TO 252 . 

300 CONTINUE 
RETURN 


f age is 
^ ^OOR qxjaljty 


SUBROUTINE OEVIA 

COMMON /AAA/X(512) , Y N ( 5 1 2 ) »YLS ( 51 2 J . OY t 5 1 2 ) , N , EM , S IGMA 
SY=0«0 - 
SQ=0, 0 

DO 100 K = 1-,N 
SY=SY+YN(K} 

DIFF=YN( Kr-^LSfK) ' ~ 

100 SQ = 5Q+(ASStDIFF) )*=5=2. 0 

SIGMA={SG/N) C.5)/ ( EY/N) 

RETURN 

END ■ ■" 


SUBROUTINE EME AN ( I NT I , INTF ) 

C■CM^ION■7AA■A7X^5l2TTy^lS^2TVYLS^BT2) Toy C512),nVEM~»S"iWaT<M ■ - - 

SYOE=OoO 

£YEDE=0.Q 

KR=0 

■28"0' 'CONTI NOE ~ - 

KR=KR+1 

K=KM-KR+1 

__IF ( K«LT.2oCR.KR.GT« INTI ) GO TO 290 

£YDE=SYdE+0'oS'*(DY(KJ+DY‘n<-l n4=(X7Ky-XlK-n ) ' ' ' " 

SYEOE=SYEDE+0«5=J*(0 Y( K ) ^-DY t K-1 ) ) -JfO » 5*=^ ( X ( K ) <-X{ K.-1 } ) * { X t K) -X ( K-l ) ) 
GO TO 280 

290 continue _ 

300 CONTINUE 
KF=KF+1 
K=KM-1-KF 

IF' <K.GT*nNVUR,KFTGT^IN‘TF)''GU TO'SIO ' " 

SYDE=SYDE + 0* DY{ K) +DY(K-1) ) * I X ( K ) -Xt K- 1 ) ) 

SYEDE=SYEDE+0.,3=5^(DY(K)+DY{K-1) ) *0 . S*( X ( X ) + X( K- 1 ) J =? ( X ( K ) -X t K- 1 M 
GO TO 300 

210 CONTINUE' ' ' 

EM=SYEDE/3YDE 

RETURN 

END 


■ SUBROUTINE S IMG ( M , N ^ EPS * lER J 

DIMENS ION A{MM} »RfNM") 

COMMON /SIM/AdOO) ,R(1G) 
NM=N*M 

■' IF'tM) 23i23, 1 ' 

1 IER=0 
FIV=0* 

DO 3 L=1,MM 

TB=ABSI'A{un' ^ 

IF (TB-PIV) 3,3,2 

2 PIV=TB 
I=L 

' 'S' "CONTI NUE" 

T0L=EPS*PI V 
l_ST=l 

DO 17 K=1,M 

IF'"{P t V >"'23723, A ■ ■ 

4 IF (lER) 7,5,7 

5 IF (PIV-TCL) 6,6,7 

6 IER=K-1 

■ 7 "FIVT=Tir7A'( I ) ■ 

J=£ I-l ) /M 

J=J+1-K 


DO a'X=K;KM,M 
LL=L+I 

TB=PIVI=^=R( LL) 

RtLL) =R(L) 

8 R'(L )=Tl3 ■ 

IF (K-M) 9.13*18 

9 L£ND=LST+M-K 
IF(J) 12*12,10 

10“TT=j!f!M ‘ ~ 

DO 11 L=LST,LEND 
TB=A(L) 

LL=I_+I I 

A(LD=ATUr) 

11 A{LL)=TB 

12 DO 13 U=LST,.Vt^,M 
LL=L+I 

’T 8"= P'l VI V5T( CCT ■ 

A(LL)=ACL) 

13 A(L)=TS 
A(LST) =J 
Piv=c'.b *’ ■ 

LST=LST+1 
J = 0 

DO 16 1I=LSTtLEND 

- prvi-^—AI IT ) 

IST=I I +M 
J = J+1 

CO 15 L=I3T.MM.M 
LL=L-J 

"A(L )=A(L T+P IVl^'ACLL) 
TB=A8S ( AtL } ) 

^ IF CTB-PIVl 15,13,14 

14 PIV=TD 

" 'I=L ■ ' " 

15 CONTINUE 

CO 16 L=KiNM,y 
LL=L+J 

16 RCLC:)*R(LL)+PIVI5!‘R(C) 

17 LST=l_ST+M 

la IF (M-1) 23,22,19 

19 IST=MM+M 

LST=M+1 • ■ 

DO 21 1=2, M 

II=LST-I 
1ST=I ST-LST 
L=IST-M 
L=A<L)+,5 

CO 21 J=II,^^/,M 
TB=R(J) 

LL=J 

DO 20 K = lST,i^f^,»v. 
LL=LL-M 

20 TB=TB-A(K5 CLL) 

— K=j+L'" : 

R ( J } =R ( K > 

21 RCK)=TB 

22 RETURN 

23 1ER=-1 

RETURN 

END 
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t 

1 


RUN PARAMETER 


DATA CARD 


READ (NSET, KMON, KDAY, KYEAR, THETAI, THETAR, PHI) 
NSET (IIP) 

> 2: Run with calibration parameters 

(or=l) (Data Card #2) 

= 2: Run without re-calibrations 

KMON, KDAY, KYEAR (3110) 

MONTH/DAY/YEAR 

THETAI, THETAR, PHI (3P10.3) 

Angular paramet ers (S.j0..»^) 

X jT 


n 

(Req^. if 
NSET > 2 
or = 1 


READ (El, M, NCVI, VI, NCV2, V2, NI, XLL, NF, XHL, INTI, 
INTF, KRG, KSTEP, NPW, NSTEP) 

El (F5 . 2) : Incidemt-beara energy 


M(I5) : Order of the lest-square fitting program 


NCVI (13) 
VI (F7.2) 
NCV2 (13) 
V2 (F7.2) 


Calibration of x-coordinate 


NI (13) 


XLL (F7.2) 
NF (13) 


Operational limits 


XHL (F7.2) 
INTI (15) 


INTF (15) 


Integration Limits 


KPG (15) ; Graphical Index 
> 1 ; with graphical output 


<1 : without graphical output 
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I 


T 

i 


n~m 

#29-if54 
ENDCAKD: 


I 

I 

i 


KSTEP (15) : Step-size for graphical output 

NPW (15) ; Output index 
> 1 : with detail output 

< 1 : without detail output 

NSTEP (15) : Step-size for detail output 

DATA. SET : For thermal background gas 

DATA SET: For reflected beam 

(Blank) 


6 ] 


